Orbitally Driven Spin Pairing in the 3D Non-Magnetic Mott Insulator BaVSs: 

Evidence from Single Crystal Studies 



O 

o 
o 

(N 

oo 
(N 



G. Mihaly^'^, I. Kezsmarki^, F. Zamborszky^, M. Miljak'^, K. Penc'''^, P. Fazekas^'^, H. Berger^ and L. Forro^ 
^ Institut de Geniue Atomique, Ecole Politechnique Federate de Lausanne, CH-1015 Lausanne, Switzerland 
^Department of Physics, Technical University of Budapest, H-1111 Budapest, Hungary 
'^Institute of Physics of the University, P.O.Box 304, Zagreb, Croatia 
Service de Physique Theorique, C.E.A. Saclay, 91191 Cif-sur Yvette Cedex, France 
^Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.B. 49, Hungary 

(February 1, 2008) 

Static electrical and magnetic properties of single crystal BaVSs were measured over the struc- 
tural {Ts = 240K), metal-insulator (Tmi ~ 69K), and suspected orbital ordering (Tx ~ 30K) tran- 
sitions. The resistivity is almost isotropic both in the metallic and insulating states. An anomaly in 
the magnetic anisotropy at Tx signals a phase transition to an ordered low-T state. The results are 
interpreted in terms of orbital ordering and spin pairing within the lowest crystal field quasi-doublet. 
The disordered insulator at Tx < T < Tmi is described as a classical liquid of non-magnetic pairs. 
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Spatial ordering of the occupancy of degenerate elec- 
tronic orbitals plays important role in the diverse mag- 
netic phenomena of transition metal compounds . To 
cite a well-known example: the interplay of magnetic and 
orbital long range ordering, and strong coupling to the 
lattice, account for the metal-insulator transitions of the 
V2O3 system In contrast, the metal-insulator tran- 
sition of the S — 1/2, electron system BaVSa is not 
associated either with magnetic long range order, or with 
any detectable static spin pairing. As an alternative, the 
possibility of an orbitally ordered ground state was dis- 
cussed [Q, while other proposals emphasized the quasi- 
one-dimensional character of the material l^-Q]. The 
crystal structure is certainly suggestive of a linear chain 
compound since along the c axis, the intrachain V-V dis- 
tance is only 2.81 A, while in the a-b plane the interchain 
separation is 6.73 A It is thus somewhat surpris- 

ing that our present studies show that electrically BaVSa 
is nearly isotropic. This means that BaVSs provides 
one of the few realizations of a Mott transition within 
the non-magnetic phase of a three-dimensional system. 
Since this case (or rather its Z? — > 00 counterpart) is 
much studied theoretically, but scarcely investigated ex- 
perimentally, a good understanding of BaVSs should be 
valuable for strong correlation physics in general. 

BaVSs has a metal-insulator transition at Tmi = 69K, 
accompanied by a sharp spike in the magnetic suscepti- 
bility [^|JlO|]. The high temperature phase is a strongly 
correlated metal with mean free path in the order of the 
lattice constant. There is no sign of a sharp Fermi- 
edge in the UPS/XPS spectra |^ and instead of a 
Pauli-susccptibility it exhibits Curie- Weiss like behavior. 
Though the magnetic susceptibility is similar to that of 
an antiferromagnet |10 no long-range magnetic or- 
der develops at the transition |^ 12|. The transition is 
clearly seen in the thermal expansion anomaly , and in 



the specific heat Q . The c?-electron entropy right above 
Tmi is estimated as ^ 0.8i?ln2, and it seems that a con- 
siderable fraction of the electronic degrees of freedom is 
frozen even at room temperature |^ . It appears that the 
69K transition is not symmetry breaking JT^ : it is a pure 
Mott transition which does not involve either magnetic 
order or any static displacement of the atoms. 

Hints of long range order were found well below Tmi, 
at Tx = 30K, in recent NMR and NQR experiments [1). 
It was suggested that orbital order may develop here, but 
it could not be decided whether the Tx phenomenon is a 
true phase transition, or purely dynamical. In any case, 
the associated entropy change must be very small since 
it escaped detection Q . In this work we prove that there 
is a phase transition, by finding its signature in static 
magnetic properties. 

In order to clarify the nature of BaVSs, we performed 
single crystal experiments, searching for macroscopic 
anisotropy in the electrical and the magnetic properties. 
Our results exclude the quasi-lD interpretations and sup- 
ply direct evidence for the dominant role ofe{t2g) orbitals 
tilted out from the chain direction. The static mag- 
netic susceptibility x, and more markedly, its anisotropy 
Xc ~ Xa, show clear anomalies both at Tmi and Tx- 

Single crystals of BaVSs were grown by Tellurium flux 
method. Powders of BaVSs and sublimated tellurium 
99.99% Ventron were mixed in a molar ratio 1 : 50 and 
heated up to 1050 °C in an evacuated sihca ampoule. 
Then it was slowly cooled down to 55°C at the rate of 
l°C/hour. The crystals, obtained from the flux by sub- 
limation, have typical dimensions of 3 x 0.5 x 0.5 mm'^. 

For the longitudinal (c-direction) resistance measure- 
ment standard four probe contact arrangement was ap- 
plied. The conduction anisotropy was determined by the 
Montgomery method [Q. Experiments performed on 
several crystals from various preparation batches showed 
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deviation only at low temperatures, which we attribute to 
the different purity of the samples (see later) . The mag- 
netic susceptibility was measured by Faraday-balance, 
the anisotropy measurements were carried out on a sen- 
sitive torque magnetometer [ p^ . 
6| . 
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FIG. 1. Temperature dependence of the resistivity Pc{T), 
and the conduction anisotropy Oc/cTa in BaVSa. The arrows 
indicate Ts, Tmi and Tx, respectively. 



Figures 1 and 2 show the temperature dependence of 
the resistivity and the conduction anisotropy of single 
crystals. BaVSs is a "bad metal" with prt — O.Tmilcm. 
With decreasing T a slight change of slope reflects the 
structural transition at Ts = 240K (Fig. 1, upper panel). 
The resistivity has a minimum at 125K. A sharp metal- 
insulator transition sets in at Tmi = 69K (Fig. 2 inset). 
Below Tmi the resistivity increases steeply and it varies 
9 orders of magnitude down to about 20K. Crossing Tx 
does not influence p{T) as it is obvious also from the Ar- 
rhenius plot (Fig. 2). We deduce a gap A « 600K for 
the insulator (A is not well defined due to a slight cur- 
vature in the In p — 1 /T plot) . The overall behavior of 
p{T) agrees well with that of high purity polycrystalline 
samples 

The conduction anisotropy, defined as the ratio of 
conductivities measured along and perpendicular to the 
chain direction, is surprisingly low, cJc/cTa ~ 3. It is tem- 
perature independent in the metallic phase and there is 
only a small jump at the metal-insulator transition. Be- 
low Tmi the anisotropy has the same small value over a 
broad temperature range down to about 30 — 40K. Note 
that in this range of T the resistivity increases about 6 
orders of magnitude in both directions. The low temper- 



ature upturn is related to impurities, the two most dif- 
ferent results obtained on different samples are plotted 
in Fig. 2 (lower panel). It seems that impurity donated 
carriers enhance the c-direction conduction, and for this 
process the activation energy is about 70K. 
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FIG. 2. Arrhenius plot of the resistivity and the conduction 
anisotropy. The latter is shown for two crystals of different 
purity. Inset: the peak of d\n p/d{l/T) defines Tmi. 

Figure 3 summarizes the results of the magnetic sus- 
ceptibility measurements. Along the chain direction, 
Xc{T) agrees well with previous data on stoichiometric 
samples; the Curie-Weiss behavior in the metallic phase 
(A'cfT = 1-2/iB; the small 9 < lOK may vary slightly with 
the range of fit) is followed by a steep decrease in the 
susceptibility below Tmi. Xc and Xa look very similar 
(therefore we show only Xc), and both are fairly smooth 
at Tx- However, we find a sharp peak in the temperature 
derivative of the a-axis susceptibility dxa/dT, and a sud- 
den break in the anisotropy Xc — Xa at Tx ~ 30K. These 
give convincing evidence of a phase transition within the 
insulating phase. We suggest that the transition temper- 
ature Tx be defined by the sharp peak in dxa/dT |16|] . 

Discussing the results first we interpret the low value of 
the conduction anisotropy. For the present purposes, we 
assume that the simple ionic picture holds (Ba^+V^+Sg^ 
ionic state, 3d^ configuration) . Figure 4a shows the crys- 
tal field splitting of the vanadium d-levels. The sulphur 
octahedra surrounding nearest neighbor vanadium ions 
are face-sharing along the c-axis. Above Ts = 240K, 
there is a trigonal distortion along the c-axis; for T < Ts, 
an additional orthorhombic component appears. The 
main effect is the trigonal splitting of the t2g level, which 
lifts the (ij,2 level (the z-axis being now the trigonal axis) 
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above the two degenerate e{t2g) orbitals |T^. The re- 
maining degeneracy is hfted by the orthorhombic distor- 
tion: below Ts, the low-energy crystal field states can be 
thought of as a quasi-doublet with a small splitting. 
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FIG. 3. Temperature dependence of the c-axis magnetic 
susceptibility Xc, the derivative of the o-axis susceptibility 
dXa/dT, and the susceptibility anisotropy, Xc — Xa- 

In Fig. 4 the lobes of the low-lying orbitals are shown 
on the background of the crystal structure. The or- 
bitals have large direct overlap along the chain direction 
as shown in the left side of Fig. 4b. Taking also into 
account that the vanadium chains are widely separated, 
any conduction through the d^^ channel is expected to be 
extremely anisotropic. In contrast, there is only weak, in- 
direct overlap of the two e{t2g) levels |0, and hopping 
along these channels is almost isotropic. From the ob- 
served small conduction anisotropy we conclude that the 
orbitals are not involved in the electron transport. 
In contrast to previous assumptions , the crystal field 
spitting between and the 6(^23) levels must be large 
and the band formed from the d^^ orbitals does not over- 
lap with the occupied states. 

We believe that the electron propagation along the c- 
axis occurs also through multiple nearest neighbor inter- 
chain hops involving only the e{t2g) orbitals. Extended 
states of the band are excited only at low temper- 
atures from impurity levels situated below the broad 
band. This process is seen in the impurity dependent 
low-temperature upturn of the anisotropy. 

The susceptibility data show that the number of po- 
larizable moments drops drastically upon entering the in- 
sulating phase. On the other hand, the metal-insulator 
transition has no magnetic precursor on the metallic side, 
X closely follows the paramagnetic formula down to Tmi. 



The volume change [|5| and the onset of an a5-plane re- 
arrangement of the V sublattice at the transition ^ in- 
dicate that the state of the system changes profoundly, 
and the effective spin-spin interaction may be completely 
different from that in the metallic phase. It seems that 
the insulator phase does not arise from a magnetic in- 
stability of the metallic phase. Let us recall the case of 
V2O3 where the onset of orbital order [|j leads to an insu- 
lating state whose magnetic ordering pattern cannot be 
anticipated from the short range order found in either of 
the neighboring phases |^ . We propose a similar picture 
for BaVSs: the metal-insulator transition involves both 
the spin and orbital degrees of freedom, and considering 
only the effect on spins, it amounts to a change of the 
spin hamiltonian. It is the consequence of the frustrated 
structure (a triangular array of V chains), and of the form 
of the relevant crystal field states, that the intermediate 
(Tx < T < Tmi) phase is not an ordered magnet, but an 
overall non-magnetic state with peculiar spin and orbital 
short range order. 




FIG. 4. (a) Crystal field splitting of the vanadium 3d lev- 
els due to the symmetry reduction from octahedral trigo- 
nal (distortion along the c-direction) orthorhombic (zig-zag 
chain developed below T5). (b) Perspective view of the crystal 
structure of BaVSa . The orientation of the low lying orbitals 
derived from t2g is also displayed. 

In order to investigate the electronic structure of the 
insulating phase, we derived a Kugel-Khomskii type 
model starting from the atomic limit, including also 
the orbital index dependence of the hopping matrix el- 
ements, and the spin-orbit coupling |l9| ]. There is a 
broad range of parameters for which the ground state 
of an isolated pair of sites is non-magnetic. Further nu- 
merical and variational calculations for clusters of up to 
24 sites revealed the presence of a large number of or- 
dered low-energy states. Figure 5 gives examples of these 
various long-period structures which obey characteristic 
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short range rules for the relative orientation of the pairs. 
Due to the small energy difference between the various 
configurations the singlet pair formation at Tmi is not 
accompanied by true long range order, but instead a 
spin-orbital liquid develops. Thermal averaging over the 
(exponentially large number of) accessible configurations 
gives rise to a homogeneous state with low susceptibility. 
Moreover, in accordance with NMR/NQR results [Q, no 
static pairing is expected over a broad temperature inter- 
val below Tmi. Long range static order is reached only at 
the much lower temperature Tx , as shown by the present 
susceptibility data, and also by anomalies in microscopic 
measurements [§,0 . Due to the pre-existing short range 
order, this transition is not accompanied by any signifi- 
cant entropy change. 
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FIG. 5. Some of the energetically favorable pair arrange- 
ments in the triangular ah plane. 

The pictures in Fig. 5 may remind us of the resonat- 
ing valence bond state of frustrated Heisenberg models 
101 . We emphasize that the situation is completely dif- 
ferent here. Within each pair, those orbitals are occu- 
pied which give rise to a strong intrapair exchange, and 
thus each cluster state belongs to a different effective 
spin Hamiltonian. Considering the shape of the e{t2g) 
orbitals one finds that pair formation does not quite sat- 
urate the exchange interaction, but there are residual 
interactions which govern farther-neighbor correlations, 
and cause a weak resonance. However, this resonance is 
much weaker than in the pure 5 = 1/2 Heisenberg mod- 
els, and the Tx < T < Tmi phase of BaVSs is better 
described as a thermal average over valence bond solids, 
than as a resonating valence bond liquid. The finding of 
a T-dependent spin gap, which vanishes at Tmi is 
consistent with our scenario. 

The metal-insulator transition affects mainly the spin 
degrees of freedom: the electronic entropy present at T 
slightly above Tmi is almost exhausted by the spin 
entropy required by the measured T > Tmi susceptibility 
which can be ascribed to ^ 70% of the V sites carrying 
nearly independent localized spins. The primary order 
parameter of the low temperature phase is the density 
of the non-magnetic pairs, and the driving force of the 
transition is the gaining of spin entropy. Sizeable orbital 
order must exist even above the metal-insulator transi- 
tion (below Tmi, it becomes additionally stabilized by the 
singlet formation as shown by the appearance of an ex- 
tra component of the orthorhombic distortion |^ ) . We 
emphasize, however, that the entropy which would cor- 
respond to the complete absence of short range orbital 



order is certainly not present even at T = 300K 0. 

In conclusion, we have shown that electron propaga- 
tion in BaVSa occurs via nearest neighbor interchain 
hops involving the e{t2g) orbitals only. The metal- 
insulator transition was described as a transition to a 
classical liquid of non-magnetic pairs, which shows spin 
and orbital short range order. We have proposed that 
though farther-neighbor inter-pair interactions prefer or- 
dered structures, due to the frustrated triangular struc- 
ture long range order can develop only well below the 
metal-insulator transition, at Tx. Our explanation is in 
overall agreement with the results of transport, magnetic 
and thermodynamic experiments. 
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